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Modified gravity has garnered interest as a backstop against dark matter and dark energy (DE). As one possi- 
ble modification, the graviton can become massive, which introduces a new scalar field. The new field can lead 
to a nontrivial equation of state (EOS) of DE which is density-and-scale-dependent. Tension between Type la 
supemovae and Planck can possibly be reduced. In voids the new scalar field dramatically alters the EOS of DE, 
gives rise to a significant gravitational slip between the two metric potentials, and also develops a topological 
defect (a domain wall) due to a nontrivial scalar field vacuum structure. 

PACS numbers: 04.50.Kd, 95.36.+X, 98.62.Sb, 98.65.Dx 



INTRODUCTION The non-detection of dark matter and 
the incomprehensibility of dark energy raises the anxiety level 
of theorists. Alternatives to Einstein's theory of general rela- 
tivity (GR) is one strategy that is useful to have should the 
cuiTent dark sector crisis persist. In modifying gravity, the 
one lesson that has been learnt is the difficulty of finding an 
alternative theory that satisfies experimental constraints on all 
scales where gravity can be probed. High precision measure- 
ments on laboratory and solar system scales must accommo- 
date astrophysical constraints on the two metric potentials in 
modified gravity theory from galactic to cosmological scales. 
Growth of density fluctuations and galaxy dynamics constrain 
the Newtonian potential. Light propagation constrains the 
sum of the two potentials. See [ 1 - 3] for review and citations. 

Modifying gravity generally introduces a scalar that com- 
plements and couples to the tensor field of general relativ- 
ity. The scalar must however be screened on small scales to 
be reconciled with the local precision tests of gravity. Three 
main mechanisms have been discussed for screening the new 
field. Chameleon f4i] and f(R) theories p"] screen via den- 
sity and decouple the scalar field at high density. Sym- 
metron fields ^ weaken the scalar coupling at high density. 
Galileon/Braneworld fields til |8|] and massive gravity (focus 
of this paper) operate via the Vainshtein effect [9J, which de- 
couples the scalar field in high-curvature regions. 

Massive gravity (by assigning the graviton a mass of order 
the inverse Hubble scale) naturally and transparently restores 
coupling on the horizon scale, where one may then hope to ac- 
count for the accelerated expansion of the universe [10]. The 
transition scale (i.e. the Vainshtein radius where decoupling 
occurs) is of order tens of Mpc today. This scale corresponds 
to the largest self-gravitating scales in the universe, and po- 
tentially the simplest to understand. 

In galaxy clusters, the metric potentials are almost the same 
as with GR, which is not necessarily tiTie in voids. The change 
in the potential is shown here to only be a few percent for 
over-densities (which has been found previously by 1 1 1-13|) 
but of order unity for voids. The apparent change in dynam- 
ics and lensing can be parameterized as an inhomogeneous 
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FIG. 1: Plot of / with Tr'/r = /(A)A^ where A = 5pni/pmO is 
the difference between the mean density inside of a radius r and the 
mean density today normalized by the mean density today. 

(scale-dependent) EOS of DE, which could be measured with 
future galaxy surveys. Even more surprising, the scalar field 
forms a topological defect (domain wall) at the edge of voids. 
MASSIVE GRAVITY Massive gravity (linearized) corre- 
sponds to a spin-2 field with five degrees of freedom, which 
can be decomposed into tensor, vector, and scalar modes. The 
effect of massive gravity in voids is well-characterized in the 
weak-field limit of GR and on sub-horizon scales, which is 
exactly the decoupling limit (see ||7l llA \lA liSl for a precise 
definition). In the decoupling limit, an effective theory char- 
acteiizes massive gravity, which only depends upon the tensor 
and scalar degrees of freedom. The vector degrees of freedom 
decouple leaving a simple Lagrangian given by lllOll . 
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The Einstein tensor {S^uhap) comes from the Einstein- 
Hilbert Lagrangian (m^j-^— g R) assessed to the second order 



in metric fluctuations /i'"^ with g^,y = ry^^ + /i^'^/mpi, the 
Newton constant SttG = nipj^, ijf^^, = (—1, 1, 1, 1), and nipi 
is the Planck mass. The second term of Eq.[T]is the standard 
coupHng of GR to the stress-energy tensor of the system. 
The third term is the Lagrangian for the scalar field tt with 

C{7T)f,. = aX^^J + {P/A^)xi^J + {-f/A^)xil\ where a,/3, 
and 7 are arbitrary parameters of the theory but are taken to 
be 0(1). Our effective scalar theory of tt is valid up to the 
cutoff scale A'^ = mpim^, where m is the graviton mass. See 
lEllstl for further details. 

Ansatz We look for solutions which can accommodate both 
under and over-densities (halos and the solar system) and are 
consistent with cosmology. The equivalence principle implies 
that locally we can always define a metric (for all time) which 
will simply be Minkowski space plus a perturbation. The per- 
turbation component will account for the non-trivial geometry 
of the space-time. Such coordinates are referred to as Fermi 
coordinates ||7l[l6|]. In the Newtonian gauge, 



ds-" 



-(1 + 2$(r))rff2 + (1 - 2*(r))da;2 



(2) 



which on small scales can be deduced from the general FRW 
metric by using the transformations tc = t+H^x'^ /2 and Xc = 
x/a{l + l/4H^x^) between co-moving (tc, Xc) and the local, 
physical Fermi coordinates {t, x). H is the Hubble parameter 
a is the scale factor We then add on local perturbations of the 
metric. $ and ^ encode both the background geometry and 
local metric perturbations. The standard stress-energy tensor 
(in Fermi coordinates) is simply T^,y ~ (p, 5ijp) where we 
have neglected off-diagonal terms (0(11^2;^) ~ u^ ~ H^z^w^ 
with velocity v). As with ^I^ and $, p and p include local 
perturbations of matter and radiation relative to a background 
density. Finally, we make the ansatz that the field tt (as well 
as the local metric and matter perturbations) are spherically 
symmetric and time-independent, tt = 7r(r). 

By the least-action principle applied to Eq[T] we obtain two 
non-trivial Equations of Motion (EOM) for the metric and one 
EOM for n with / = d/dr 
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aV2(2vI/-$) + M^V2(vI/-$)-5i^V2<I' = (5) 
^ ' A^ r ^ ' A^ r^ 

where in Eq. |4]we are not summing over i but indexing over 
the spatial dimensions {x,y,z). Previous authors lll2i 11311 
have found the above EOM, but have neglected the effect 
of pressure in Eq. |4] Also, ( ) gives the average value of 
a given quantity inside of a radius r. The quantities (p) 
and (p) are the average values of p and p inside of a radius 
r (with contributions from matter, DE, radiation, etc.) and 



(p,^) = M(r)/(4/37rr3) where M{r) = j p^{r')r''^dr'dn 
and Pm(?') is the density only in matter at a point r. 

In linearized GR, perturbations and background evolution 
can be separated out by simply subtracting off the background 
evolution from the perturbed Einstein equations. Naively one 
would expect the same in massive gravity, but this proves im- 
possible. Eq. |5] is non-linear in tt and in the metric pertur- 
bations, which introduces cross-terms between the perturbed 
solution and the background solution. At zeroth order in tt we 
look for solutions that accommodate both the background and 
local perturbations in the metric, density, and pressure. 

The above EOM can be made dimensionless by multiplying 
by m^^ with A"^ = mp\m^ where m is the graviton mass and 
then by setting tt' /r = /(A) A'^ where A ~ Spm/pmo with 
p„iO ^s the matter density today and Jp^ is the under-or-over- 
density of a void or a halo with Jpm = (pm) ^ Pmo where 
(5p,n and (pi„) depend upon r. Hence, 5p„x/ PmO — ^^ cor- 
responds to an empty void as in devoid of matter but not DE. 
After making the appropriate substitutions, Eq. |5]becomes a 
quintic constraint equation for /(A) which depends upon the 
average pressure (p) and density (p) for a fixed radius r. Out 
of the five roots, three will typically be real. 

Classes of Solutions The solutions for tt' can be catego- 
rized into three separate classes. In a separate publication we 
will discuss in detail the various solutions, and only outline 
the different solutions here. Two of the classes have been dis- 



cussed previously by Ill2lll3ll . Neither of these cases produce 
interesting cosmological solutions. One class of solutions de- 
gravitates all mass. The second set of solutions generates an 
EOS of the universe which is equivalent to radiation at late 
times and generates negative densities. 

The third class gives a new cosmological solution in which 
the EOS of DE generally tracks the energy density of matter, 
radiation, and curvature. There are several different sub-cases. 
In the first self-accelerating case (no cosmological constant), 
we can tune a, /3, and 7 to effectively generate a cosmological 
constant. The authors of 1 10] found a similar solution but with 
a different ansatz, which decoupled the evolution of the scalar 
field from the evolution of matter, radiation, and the metric. 
In a second self-accelerating solution (in a flat background), 
the EOS of DE will depend on the local density of matter 
In a third self-accelerating case (with curvature k ^ 0), the 
scalar field can dynamically counteract the effect of curvature 
upon the expansion rate of the universe, which may appear as 
so-called "phantom" EOS of DE with w < —1. 

The more general case with a cosmological constant will 
constitute the focus of the paper In solving for the tt field 
EOM, we search for roots of a quintic equation which has no 
general analytic solution. We search for all solutions numer- 
ically, which requires picking a, /?, 7 and m (graviton mass) 
such that densities are positive, the field tt is ghost-free, and 
constraints on Hq and the f^m are satisfied. We will discuss 
the overall parameterization in a separate publication, but to 
be concise we will consider a single case with a — —0.45, 
/3 = —1.0, and 7 = —0.86 in which a topological feature ap- 
pears (for other choices of parameters no domain wall exists). 



The same solution exists for fixed ratios of 7^/a and 0^ ja. 
Finally to satisfy cosmological constraints, we have also set 
^irGpfj/'m? — 7.65 where pg is the measured critical density 
today and have set rj,n = 0.3 (the qualitative results presented 
are invariant under varying fim and Hq). 

In Fig. 1, we have plotted the solution of vr'/rA'^ = /(A) as 
a function of the average under- or over-density 5pai interior 
to the radius r over the average density in matter today pm 
with A = (^Pm/Pm ■ Fig- 1 shows how the field changes from 
the case of an under-density to an over-density. 

Over-Density First, we consider a point mass A/ > (as 
our over-density). In the limit of approaching the point mass 
r ^ 0, then A ^ 00 and /(A) -^ -0.42. The Vainshtein 
mechanism kicks in recovering GR($ — ^ 'J -> —GM/r). As 
r ^f CO then we smoothly go onto the cosmological solution 
and /(A) takes on a value which gives the correct cosmologi- 
cal solution today with 2$ -H- -(H + W^)x^ and 2* -^ H^f^ 
where H is the Hubble parameter. From the metric EOM, we 
then also recover the Friedmann equations. 

Under-Density Far outside the void 5p^ ~ (the metric 
potentials match onto the cosmological solution found in the 
previous section). The average density (pni) drops as we move 
to the center of the void. Eventually at a certain radius, the 
derivative of the field must jump to a new branch (see Fig. 1) 
once (5/9in drops to less than half of the present mass density 
of the universe today. The jump forces a discontinuity in the 
derivative of the field tt' i.e. a kink, but in principle, we have 
the freedom to adjust the integration constant of vr' such that tt 
is continuous. Regardless, to match the derivative across the 
discontinuity requires the introduction of a wall. 

The nature of the tt wall is categorically different from a 
wall formed for a typical scalar field. A domain wall typically 
corresponds to a kink solution interpolating between two dif- 
ferent vacuum states on different sides of the wall. The tension 
in the wall is due to the potential separating the two vacuum 
states in field space. The tt field has no potential, V{tt) ~ 0, 
but upon removing the static condition, the EOM for tt be- 
comes highly non-linear As a parallel in fluid mechanics, 
non-linearities of the Navier-Stokes equation can lead to a dis- 
continuity in density and pressure gradients (which we inter- 
pret as a shock wave). Similarly for vr, a "shock wave" could 
also dynamically generate a discontinuity in vr'. The presence 
of vacuum energy and matter are critical to induce the appear- 
ance of the domain wall. The vr field in empty space will not 
induce topological structure (see lll7ll for a proof). 

From Fig. 1, there exists a mapping of tt' /r into real space, 
which effectively maps the vacuum of vr into real space, 
since there is a one-to-one correspondence between vr and vr'. 
Clearly the mapping of tt' is not one-to-one into real space, but 
we speculate that in a higher-dimensional space such a prob- 
lem might be avoided. We can amend this problem by noting 
that a quintic equation corresponds to a curve on a Riemann 
surface which has multiple branch points. In fact, / loops 
around two branch points producing the multivalued function 
in real space in Fig. 1. On the Riemann surface, / is single 
valued! There exists a one-to-one mapping between real space 
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FIG. 2: Plot of EOS of DE iude where {pm) / p-mo is the mean 
density of matter interior to a radius r over the mean density today 

PmO- 

and 7r'/r for / on the Riemann surface. If we allow a complex 
extension of A and vr', the quintic equation corresponds to 
hyper-elliptic surfaces which are compact Riemann surfaces 
with a genus greater than one which we identify as the vac- 
uum of vr by analogy with a standard Nambu-Goldstone Bo- 
son (NGB). 

By the above identification, we can make sense of the ap- 
pearance of a domain wall. The vr field can be interpreted 
as a NGB (see fllSll '). Frequently the topologically non-trivial 
vacuum of a NGB when mapped into real space leads to the 
appearance of topological structures. Similarly, our function 
/ due to its wrapping around branch points is equivalent to 
a non-trivial wrapping around the hyper-elliptic surface. Fur- 
ther properties of the domain wall may potentially be inferred 
from the Riemann surface. 

Dynamical Dark Energy Massive gravity effects DE. The 
mean density inside a sphere of radius r is (/3(r)) — {p^{r)) + 
Pec where pm is the cosmological density of matter today, pcc 
is the energy density from the cosmological constant. Massive 
gravity also contributes to the energy density of the universe 
(Pmg) via the last terms of Eq. [3] Similarly in Eq. |4] (p) — 
— Pec neglecting radiation. Massive gravity also adds to the 
pressure of the universe p^g via the last terms in Eq. |4] 

Massive gravity naturally leads to inhomogeneous DE 
evolving in time. In Fig. 1, the vr field responds to changes 
in the local density. As vr evolves, so does pmg and p,„g. The 



EOS of DE is wde ~ 



{P)+Pn 



We have plotted wue in Fig. 2 



as a function of the average density over the mean density to- 
day. Today (outside of halos and voids) wue = —0.97. As we 
move into a halo, wde asymptotes towards —1. Conversely, 
at the center of a very under-dense void, wde jumps to —0.69. 
We can also read Fig. 2 as a measure of how the EOS of DE 
changes with redshift. At high redshift, the average density 
of the universe becomes large compared to the mean density 
today, equivalent to taking (pm)/pmO — ^ oo- At early times, 
the scalar field becomes unimportant which is merely a man- 



ifestation of the Vainshtein mechanism. In this epoch, DE is 
then unimportant for the evolution of the universe. 

Recent cosmic microwave background results from the 
Planck satellite [ 19] indicate (with 2.5a significance) a nearly 
10% lower value of the Hubble constant (Hq) compared to 
values found from local-Universe and supernova-based mea- 
surements, e.g. [20'|. Although a number of SNIa systematics 
and concordance-model uncertainties could account for this 
difference, such an effect could also be explained within our 
model by a significant fraction of low-density void in the local 
Universe or along the line of sight. For a difference Aw in the 
background EOS, AH/Hq « SzAw/2. With a difference Aw 
driven by wde ~ —0.7 in low-density voids and wde ~ — 1 
elsewhere, a higher value of Hq would be inferred than from 
the cosmic microwave background. The required local den- 
sity of < /Oi no/ 2 ap pears consistent with recent observational 
studies, e.g. fl2Ul22ll . This could potentially relieve the tension 
between the different Hubble constant estimates up to the ob- 
served 10% level. 

CONCLUSION Massive gravity can naturally give rise to 
an inhomogenous EOS of DE, which will be dependent upon 
local matter density and scale. The scale-dependent EOS of 
dark energy may reduce the tension between relatively local 
measurements of dark energy from SNIa and global measure- 
ments from the CMB. Detailed analysis of voids may be the 
most direct way to infer the existence of a scale-dependent 
and density-dependent EOS of DE since perturbations in mat- 
ter will lead to perturbations in DE. We can infer the exis- 
tence of inhomogeneous DE via dynamics and lensing. The 
motion of galaxies from the geodesic equation constrains $ 
(in the weak-field limit). However, for massless particles, the 
geodesic equation (lensing) depends upon the lensing poten- 
tial $L = l/2($ + 4'). We can infer the amount of material in 
and around a void by simply counting galaxies. It is then pos- 
sible to infer the EOS of DE in a void and test the predicted 
large gravitational sUp between the metric potentials, which 
should be achievable in future surveys 1231 12411 . 

The EOS of DE will change the growth rate of voids. The 
ISW effect can also be a good way to constrain the EOS of 
DE in a void |25]. The domain wall may also directly af- 
fect lensing, dynamics of galaxies and cosmic rays crossing 
the domain wall. With a magnetic field, charged particles 
could be accelerated via an interaction with the tt field, and as 
with axions/pseudo-Nambu-Goldstone bosons, distinct pho- 
ton conversion/polarization could take place. 

Dimensionally, the tension in tt goes like t ^ A^. Match- 
ing the metric potentials across the the discontinuity in tt' will 
give a precise solution for the tension in the domain wall. We 
should also be able to calculate the wall tension via non-linear 
field dynamics and via an analysis of the vacuum structure of 
TT. GR implies that they should all give the same solution 



renormalization theorem may be a consequence of the com- 
pact vacuum of tt, as was found in the original WZW model. 
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Finally, the non-renormalization theorem for galileons 112 611 
appears mysterious. The authors of 1181 showed that the 
terms in C{tt) are in fact Wess-Zumino-Witten (WZW) 
terms. If indeed our solutions for tt are 'onto', then the non- 
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